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Abstract 

We study numerically on the lattice the 2D Yukawa model with the U(l) chiral symmetry 
and Np = 16 at infinite scalar field self-coupling. The scaling behaviour of the fermion mass, 
as the Yukawa coupling approaches zero, is analysed using the mean field method. It is found 
to agree with that of the Gross- Neveu model with the same symmetry and Np. This is so 
even if the sign of the bare kinetic term of the scalar field is negative. This suggests that 
the 2D Yukawa models belong to the universality class of the Gross-Neveu models not only 
at weak scalar field self-coupling but also for a broad range of the bare parameters which is 
not accessible to the 1/Np expansion. New universality classes might arise at the crossover 
to the spin model universality class, however. 
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1 Introduction 



The 2D Gross-Neveu (GN2) models with various global chiral symmetries [|l[] have been of 
continuous interest since it has been realized that these models are asymptotically free and 
that the dynamical mass generation (DMG) occurs ||I|, [J. The 2D Yukawa (Y2) models with 
the same symmetries have received very little attention, however, though they can potentially 
have very similar properties. These models are defined in the continuous Euclidean space 
by the action 



2 ^ 4! 



where mo, A and y are the bare scalar mass and bare dimensionful coupling constants, 
respectively. For definiteness we have written down the model with chiral U(l) symmetry 
and a quartic scalar field selfcoupling. The generalization to other symmetries and/or other 
selfcouplings is obvious. The fermion field ip consists of Np flavour components and if is a 
complex scalar field. 

It was noticed long ago that the Yukawa models in various dimensions reduce to 
the four-fermion models if the bare kinetic and selfinteraction terms for the scalar field 
vanish. This relationship in 4d was used in the models of composite Higgs boson M. Such 
a relationship is a rather simple kinematical fact which is even more clearly seen in lattice 
regularized Yukawa theories |5|, ||] . 

In this paper we present numerical results suggesting in 2D a much more profound re- 
lationship, namely that the Y 2 field theory belongs in a broad range of parameters to one 
universality class, which is the class of the GN 2 model with the same symmetry. Such equiva- 
lence might be indicated by the dimensionality of the coupling parameters in the Y 2 models, 
but we are aware of no clear demonstration on such a general level. For large Np and small 
y, A = 0(1/Np) the equivalence can be shown in the large Np expansion Already fifteen 
years ago the equivalence between Y 2 and GN 2 has been suggested by Guralnik and Tam- 
vakis m on the basis of a mean field (MF) analysis for a wide class of scalar selfcouplings. 
However, the MF method is traditionally considered not to be reliable in 2D. 

Recently we have found numerical evidence that in analogy to the GN 2 models the Y 2 
models with the Z(2) and U(l) symmetry on the lattice are asymptotically free in the Yukawa 
coupling even if the selfcoupling of the scalar field is arbitrarily strong || |K]] . This property 
is suggested by the MF method [K| and in the effective potential approach |12| which is 



in spirit similar to the MF approximation. It is argued in ref. ||T0[ that the MF method 
is applicable to the Y 2 models because the fermions induce long range effective interaction 
which couples many scalar variables ferromagnetically. Such interaction can be well described 
by a mean field. The local selfinteractions of the scalar field (without fermions) are treated 
exactly or numerically. Such a version of the MF method is applicable on finite lattices 
and gives a good agreement with the numerical data. This allowed us to demonstrate the 
exponential decrease of the fermion mass with 1/y 2 , and thus the asymptotic freedom of the 
Y 2 models. 

The agreement between the numerical and MF methods obtained until now suggests to 
test numerically for various values of the coupling parameters and with improved data the 



MF inspired conjecture [@, [L2|, [HJ that the Y 2 and the GN 2 models with the same symmetry 
belong to the same universality class, i.e. correspond to the same renormalized quantum 
field theory and have the same physical content in a broad range of the coupling parameters. 

For this purpose we have investigated quantitatively the rate with which the fermion 
mass amp in lattice units (a being the lattice constant) approaches zero as the bare Yukawa 



coupling y decreases. We compare this ^-dependence with the MF predictions for various 
fixed values of the dimensionless bare parameters a 2 m^ and a 2 A . The parameters should be 
fixed such that the pure scalar model at y = is in the high temperature phase (symmetric 
phase in the Z(2) case, vortex phase in the U(l) case, etc.). The expected asymptotic scaling 
law in the infinite volume is 

am^-exp^-Li) , (1.2) 

where 

g 2 = xv\ (1.3) 

y is the bare dimensionless Yukawa coupling of the lattice Y 2 model and \ is the susceptibility 
of the pure scalar model at y — 0. For positive values of the hopping parameter k of the 
scalar field (its definition is given below) one has 

x = (1.4) 

y 2 = 2na 2 f, (1.5) 

am^ being the renormalized mass of the scalar boson in the high temperature phase and 
its wave function renormalization constant. 

We find numerically that in the U(l) model with Np = 16 the coefficient /3q is within 
certain accuracy independent both of am^ and of the strength of the bare selfcoupling of the 
scalar field. This is so even if this coupling approaches infinity. The values of (3q obtained 
from the numerical simulation are close to the value of the first coefficient of the /3-function 
of the corresponding GN 2 model. The agreement is good (within 7%) at large am^ and 
is found to be even better if the sign of the bare kinetic term of the scalar field is chosen 
negative. On the other hand, for small am^ significant deviations from the MF prediction 
are observed in the range of couplings and lattice sizes we have investigated. This might 
indicate a crossover towards the scaling properties of the pure scalar theory. 

The parameter region corresponding to the negative sign of the bare scalar kinetic term 
(negative k) is of particular interest. Here, far outside of the scope of the perturbation 
expansion, the contributions from the fermion determinant to the effective kinetic term of 
the scalar field make the Y 2 models apparently field theoretically as consistent as the GN 2 
models. 

Further support for the conjecture of the universality of the Y 2 and GN 2 theories is 
provided by our numerical results for the ratio of the mass of the scalar fermion-antifermion 
bound state B and of the fermion mass, tub/wif- The mass am# is extracted from the 
two-point function of the composite operator ipip. The ratio is found to be independent of 
the chosen values of the bare parameters of the Y 2 model. For Np = 16 fermions it is close 
to 2 and in agreement with the predicted [O, KM value 1.99 in the GN 2 model. 



The outline of the paper is as follows: In the next section we define the lattice Y 2 models 
and describe the structure of their phase diagram, as implied by the MF method and our 
previous numerical results. The various universality classes, suggested by this structure, are 
discussed. In sec. 3 we show how the universality of the Y 2 models can be studied by means 
of the scaling behavior of the fermion mass. Here we also describe how this scaling behavior 
can be determined, using the MF method, on lattices of modest sizes and at rather large 
values of the Yukawa coupling. Section 4 presents numerical evidence obtained within the 
U(l) symmetric Y 2 model with Np = 16 that the scaling behavior of the fermion mass is very 
close to that of the GN 2 models. Some differences are noted too, however, and interpreted as 
a sign of a cross-over to the spin-model universality classes. In sec. 5 we further support the 



idea of equivalence between the Y2 and GN2 models by comparing the ratio of the two lowest 
masses. Section 6 and the appendix contains some analytic speculations about the cross-over 
between the GN 2 and spin-model universality classes. A brief summary is presented in sec. 7. 



2 Phase diagram of the Y2 model 

We investigate Yukawa models with Z(2) and U(l) chiral symmetry using mostly staggered 
fermions. The lattice actions have been chosen such that the GN 2 and the spin models are 
their simple special cases. E.g. the action of the U(l) model for staggered fermions is: 



S = S B + S F + S Y (2.1) 



i-i- 



(2.2) 



x,a,fj, 

= y E5£ E (<f>l* + i(-l) Xl+X2 <t>l- b ) X a x (2-4) 

x,a beP 

where <fi is the two-component scalar field and x° S a = 1, - - - , A r , are N staggered fermion 
fields, which describe in the continuum limit Np = 2N Dirac fermions. We chose the 
representation rji x = (— 1) X2 and r]2 X = 1 for the sign factors. The dimensionless bare 
parameters are k, the scalar hopping parameter, A, the scalar quartic selfcoupling and y, 
the hypercubic Yukawa coupling, which couples the scalar fields on a plaquette P to the 
fermionic fields located in one corner of P. The action for the Z(2) model is similar, cf) are 
then only 1 -component scalar fields. 



The continuum action (|1.1|) arises from the lattice action ( |2.1|) if a — > and the following 
rescalings are made: 

ml = (1 - 2A -4/«)/(a 2 /t) 
A = 6A/(aV) 
y = y/(aV2n) 



(2.5) 



The components of ip arise from those of \j y/a in a complicated but standard way. 

In the case k = A = this action describes the Gross-Neveu models on the lattice with 
4> being an "auxiliary" scalar field. The usual GN 2 coupling constant is 

9 = y/V2. (2.6) 

The parametrization used in ( |2.1| ) makes it obvious that the GN 2 models are special cases 
of the Y 2 models. In the parametrisation usual in the continuum ( |1.1|) the case k = A = 
corresponds to the nonperturbative limit a 2 ml, a 2 y 2 — > oo, y 2 /ml = g 2 fixed. 

The other well understood limit is y = 0, where the action ( |2.1|) describes the pure 2- 
component scalar field theory and free massless fermions. In the limit A = oo it is the 2D 
XY model. In the Z(2) case the 2D Ising model is obtained in this limit. 

The 1/Nf expansion usually used for the GN 2 model can be applied to the Y 2 models 
only for A = 0(1/ Np) and leads to results which are equivalent, at least in the leading 1/Np 
order, to the GN 2 case 0. In [0 we argued that the MF approximation can be used to 



describe the scaling behavior and also some finite cutoff and finite size effects in the Y2 
models for any A. The fermion determinant generates nonlocal ferromagnetic interactions 
between 0's. The MF approximation is applicable because a large number of sites contribute 
to the mean field. 

The method is briefly as follows. The effect of the fermion determinant can be described 
by a mean field H(a, y) acting on as an external field, 

H(a,y) = Ny'a J ^ fesin 2 ^ + (yaf\ , (2.7) 



where a is the mean magnetization. The mean field ( |2.7|) acts as an external field on the 
scalar model given by the action Sb from ( |2.1| ) and the magnetization (0) is given by the 
selfconsistency equation 

(0) = /(if«0),y)), (2.8) 

and amp = y(<p). Here f(H) is the response function of the scalar model. On a finite 
lattice the momentum integrals in H((<j)),y) are replaced by sums over lattice momenta. 
The resulting equation is solved numerically by recursion predicting (0) as a function of the 
Yukawa coupling y, the scalar couplings k, A and the volume of the system L 2 . In order 
to obtain the response function f(H) we study the scalar model, with interactions from the 
pure scalar part of the action, in the external field H. This we can do numerically and in 
some special cases exactly. The mean field H can then be determined selfconsistently. 

We note that the logarithmically diverging nonlocality of the ferromagnetic interaction 
producing the mean field in ( |2.7| ) makes the DMG possible at arbitrarily weak Yukawa 
coupling and leads to the asymptotic freedom. 

As discussed in [IT] the MF approximation leads also to a better understanding of the 



fermion mass generation mechanism in the case of the Y 2 model with continuous chiral sym- 
metry which cannot be broken spontaneously according to the Mermin-Wagner-Coleman 
theorem. The main idea is that the effective interaction induced by the fermions can produce 
ferromagnetic ordering of scalar variables on distances 0(l/mp). On larger distances the 
ordering is destroyed by the long wavelength fluctuations (spin waves) as usual in two dimen- 
sions. The scalar field in a volume of a linear size 0(1/ mp) has a nonzero magnetization 
(4>)i/m F which direction is drifting slowly. The mass of the fermion coupled to such a scalar 
field is given by am F ~ y(<P)i/ mF - We checked this relation and found a good agreement in 
our Monte Carlo simulations. One can view this relation as a selfconsistency equation for 
mp. Alternatively it can be used as a formal definition of (0)i/ mF . Thus we can apply the 
MF approach to the Z(2) and U(l) models in complete analogy if we bear in mind that by 
a we mean (0)i/ mF . 

Figure [I] is a schematic representation of the phase diagram of the Y2 models with Z(2) 
or U(l) chiral symmetry. At y = we have a critical line k c (X) separating two phases of 
the pure scalar field theory. In the Z(2) case these are the high temperature, paramagnetic 
(PM) phase and the low temperature, ferromagnetic (FM) phase. In the U(l) case the 
corresponding phases are called vortex (VX) phase and spin wave (SW) phase, respectively, 
alluding to the dominant spin configurations. 

There is considerable evidence provided by some exact, semiclassical and large Np results 
[p], 13, |14], H|, [H| that the GN2 models are asymptotically free and exhibit the DMG. This 



implies that on the line k = A = the systems are for arbitrarily small coupling y in the 
broken (FM) phase for the Z(2) symmetry or in the SW phase in the U(l) case. The models 
have an essential singularity at y — 0. Analytic and numerical investigations of the Yukawa 
models || || [12|, |l0j indicate that for n < k c (\) the theories behave for any A similarly to 




Figure 1: Sketch of the phase structure of 2 dimensional Yukawa models with hypercubic 
Yukawa coupling in the 3 dimensional coupling parameter space. 



the GN2 models. They are asymptotically free. At arbitrarily small y the models are in the 
FM or SW phase for the Z(2) or U(l) case, respectively. The surface k < « C (A), y = is a 
critical surface at which the masses of the GN 2 model vanish, if expressed in lattice units. 
One can perform the continuum limit approaching any point of this surface. For some lattice 
discretizations of the GN 2 or Yukawa models on the PM or VX phases can appear at large 
y (16], [J. This is a lattice artifact, which can be avoided by using the hypercubic Yukawa 
coupling. 

The scalar part of the action ( j2.1|) has the so-called staggered symmetry 

<P X - (-ir +X2 x =: k - k. (2.9) 

In the scalar models this symmetry implies a phase transition line at K — — k c (X) which 
separates the PM and the antiferromagnetic phase in the Z(2) model and the VX phase 
from the staggered spin wave phase in the U(l) case. 

The staggered scalar field <p st does not couple to the fermions through the hypercubic 
Yukawa coupling ( |2.4| ). Therefore we expect the phase transition at k — —k c (X) to occur 
also for y > and to depend only weakly on y. It forms a critical sheet which separates in 
the Z(2) case the FM from the FI (ferrimagnetic) phase (see fig. |l|). We have checked the 
existence of the FI phase in the case of the Z(2) model numerically. A similar critical sheet 
is expected also in the U(l) model. 

Though we do not know whether the critical sheet y = 0, k < k c (X) extends to k — —00 
or not, we see no reason why it shouldn't. The phase transition around k ~ —k c (X) does 
not influence amp which seems to scale with y — > also below this sheet. 

The described phase structure suggests the existence of several universality classes in 
the Y2 models. They depend on the position of the critical point where one performs the 
continuum limit and, as we shall see, in some cases probably also on the way how this point 
is approached. 

When one performs the continuum limit on the critical line k c (A) (A > 0) and approaches 
it within the plane y — 0, in the Z(2) case the theory belongs to the Ising universality class. In 
the U(l) case the transition is Kosterlitz-Thouless-like, the theory belongs to the universality 
class of the 2D XY model. At A = and k = 1/4 there is the Gaussian fixed point. When 
approaching this point from smaller k at y — we have a theory with free noninteracting 
scalars and free massless fermions. The continuum limit towards the line k = k c (\) taken 
from the y > region is rather complex and will be discussed in section 6. 



Approaching the critical surface y — 0, k < k c (X), leads to an asymptotically free theory 
H ^, |12|, But it is not a priori clear whether this theory is for any A > and n < k c (X) 
in the same universality class as the GN 2 model at A = k = 0. At least the spectrum is 
expected to be similar, with fermions and fermionic bound states scaling while the masses 
of the states in the 4 theories, in particular am^, remain nonzero and decouple in such 
a continuum limit. The discussion of this limit is the main aim of the next two sections 
of this paper. We shall concentrate on two aspects of universality: the question whether 
the theories scale in the same way when one approaches the critical surface at different 
(k, A) points. More precisely, whether the first coefficient of the /3 function, which can be 
interpreted as a critical exponent, is the same. We also examine the ratio of the first two 
masses in the spectrum of the Y 2 model and compare it with that of the GN 2 model. 

We note that the structure of the phase diagram and the continuation of the critical 
surface at y = to negative k suggests to consider also k < 0, though the rescalings ( |2.5|) are 
not well defined. Even if negative k corresponds to the bare scalar kinetic term with wrong 
sign, the continuum Y 2 theories can be physically sensible at k < 0. A strong indication for 
such a phenomenon has been found numerically in the 4D Yukawa theories on the lattice 
for sufficiently strong Yukawa coupling |T7|]. Thus we shall study the scaling behavior at 
negative k also. 

At the critical surface around /t ~ — /t c (A), y > 0, the fermion mass does not scale and 
fermions would decouple in the corresponding continuum limit. 



3 Scaling behaviour of amp and its analysis on finite 
lattices 

The asymptotic scaling law for the GN model (1 loop) is 

am F ~ exp ■ (3-1) 
We note that the correct values of /3q obtained in perturbation theory are 

Po = ^-± (Z(2)), (3 = ^ (U(l)). (3.2) 
Zix Lix 



The MF considerations [ |TU| , 15 ] and the study of the effective potential |T2| predict that 
the fermion masses in the more general case of the Y 2 models behave according to the scaling 
law (p..2[). Here x — x( K > is the susceptibility of the pure scalar model described by the 
part Sb of the action ( |2.1| ). For < k < k c (A) it is given by (|L4]), but the scaling law ( |1.2| ) 
holds also for k < 0. Indeed, x is finite for such k values. In particular, 

x(k,0) = 1 —^. (3.3) 
Ik 

The predicted value (3q = N f /2ti does not depend on k and A and thus coincides with the 
value of the first coefficient of the (3 function of the GN 2 model in eq. (|3.1| ), as predicted 
by the MF approximation in that case. In the U(l) case the MF approximation gives /?o 
correctly. We note that the scaling law ( |3.1|) of the GN 2 model agrees with (|1.2|) because of 
eq. (p|) and X (0,0) = l/2. 

As mentioned in the previous section we concentrate on the question whether the con- 
tinuum limit taken at the critical surface y = 0, k < k c (X) is everywhere Gross-Neveu like. 
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Figure 2: The 4 point function in a Yukawa theory in the lowest order in y and the scalar 
selfinteraction fully included. 

Therefore we want to compare the scaling behaviour of the fermion mass amp for y — > at 
various fixed k and A. In order to make such a comparison it is necessary first to compare 
the strength of the Yukawa coupling itself. The interacting scalar field namely modifies 
the effective strength of the Yukawa interaction already on the fermionic tree level. This 
is obvious when the fermion four-point function Gp is considered in the lowest order in y, 
namely 0(y 2 ), but with the scalar selfinteraction fully included, as shown in Fig. |^. At zero 
momentum we obtain 

G?( Pl = 0)=y 2 Gf(p = 0)=y 2 X . (3.4) 

Thus we see that the scaling behaviour should be compared with respect to the effective 
strength of the Yukawa coupling 

g 2 = xy 2 - (3.5) 

(We thank F. Niedermayer for an elucidating discussion on this point.) Therefore the scaling 
behavior of amp for various fixed k, A values is interpreted as the same if the numerically 
determined coefficient (3$ in the scaling law ( |1 . 2|) is the same. Its consistency with the values 
( |3.2| ) is then an indication that the Y 2 models belong for those k, A values to the GN 2 
universality class. 

The aim of our numerical investigations is thus to check whether the scaling law (1.2) 



is consistent with the data. However, eq. ( |1.2j ) is valid only in the infinite volume limit at 
very small y. In ref. [[H|] we described how the MF method can make predictions about the 



fermion masses on finite lattices. 

These MF predictions allow us to estimate how Yukawa models approach the asymptotic 
scaling. In fig. |3| we have plotted the fermion mass amp obtained by solving eq. ( p.8[ ) 
for k, = A = on lattices with L = 16, 32, 64 and 256. The onset of asymptotic scaling 
represented by the straight line is around 1/g 2 ~ 20, i.e. it can be seen only on very large 
lattices (L>200). 

Figure § shows the MF prediction for amp at A = oo and k = —0.2. The picture is 
similar to the GN 2 case, but an indication of the same asymptotic scaling seems to require 
even larger lattices. For other k at A = oo the situation is analogous. 

The MF analysis thus suggests such a slow approach to the asymptotic scaling that we 
have no chance to achieve it in numerical simulations. Too small y and, correspondingly, 
too large lattices would be needed. The strategy we adopt in this situation is the following: 
we compare the numerical data for amp with the MF predictions for finite volumes and 
correspondingly large y. If we find an agreement within a reasonable accuracy margin, we 
conclude that the MF method works well and that its asymptotic predictions are thus correct. 
In other words, the MF analysis of the data is the theoretical means for extrapolation to 
large 1/g 2 and L of the data obtained numerically at moderate values of these parameters. 

To compare numerical data with the MF predictions we introduced an ad-hoc correction 
r to be determined in fits to the data. It's role is to take into account eventual deviations 
from the mean field approximation and to quantify the (dis) agreement between data and 
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Figure 3: Mean field approximation prediction for the fermion mass in the GN case, plotted 
logarithmically against 1/g 2 = l/(xy 2 )- The slope of the straight line corresponds to the 
asymptotic scaling, eq. fl3.1|) . 
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Figure 4: Same as fig. |, now for A = oo, k = —0.2. 
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Figure 5: Numerical data for fermion mass nip on finite lattices in the U(l) Y 2 model plotted 
against 1/g 2 , where g is the effective Yukawa coupling g = ^fxy. 



the prediction. It is introduced in the following way: instead of amp = y{4>) we use in eqs. 
( |2.7| ), ( |2.8| ) amp = ry(<f)) and fit for the parameter r. In this MF motivated Ansatz the 
asymptotic scaling would have a different form, namely 



amp ~ e ^ca 7 . 



However, on lattices much smaller than the sizes at which the asymptotic behavior sets on, 
r parametrizes the difference between the MF prediction of the approach to the asymptotic 
scaling and the true one. So we expect r — > 1 as L — > oo and y — > 0. The result r = 1 would 
mean perfect agreement between the data and the MF prediction for the coefficient (3$ in 



the asymptotic scaling behavior in eq. ( |1.2| ), which of course cannot be expected on finite 
lattices. 



4 Observed y dependence of amp 

In order to verify the hypothesis that GN and Y 2 models are in the same universality class 
as y — > for k < k c (A) we have studied the ^/-dependence of am F in the U(l) model with 
Np = 16. The numerical data have been collected using staggered fermions on lattices of 
the sizes L = 16, 32, 64 at many y points for five pairs of fixed k and A values (see table 1). 
Using the hybrid Monte Carlo algorithm, the fermion propagators have been measured in 
coordinate and in momentum space and have been fitted with the free fermion Ansatz. The 
results for amp are very stable and the errors so small that the ^-dependence of this mass 
clearly shows deviations from the MF predictions on finite lattices. 

The pair n = A = (the GN 2 model, which is the best understood case) has been 
chosen to estimate the quality of the MF description of the data. For this purpose we have 
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Figure 8: Same as fig. [j[ for A = oo and k = 0.2. 



k A 


16 2 32 2 64 2 



-0.4 oo 
-0.2 oo 

oo 
0.2 oo 


0.957(3) 0.956(5) 
0.971(5) 0.970(7) 
0.955(5) 0.95(1) 0.94(3) 
0.94(1) 0.93(2) 0.93(4) 
0.91(4) 0.86(4) 0.85(6) 



Table 1: Values of the fit parameter r in the U(l) Y2 model with Np = 16 on 16 2 , 32 2 and 
64 2 lattices. 



determined the value of the parameter r by a fit to the data as described at the end of the 
preceding section. The fit is shown in fig. 5 and we find that r has the same values for both 
lattice sizes L = 16, 32 we have used in this case and deviates from 1 by 4.5%. Having chosen 
rather large N F = 16 we can be sure that the GN 2 model scales according to eq. (|1.2j ), as 
predicted by the large Np expansion. So the deviations of r from 1 of this magnitude should 
not be interpreted as a signal for a scaling behavior different from (1.2), but rather as an 
estimate of the accuracy within which the MF method works on finite lattices. 

In the figs. H and |8] we present numerical results of the hybrid MC simulation of the Y 2 
model for A = 00 and k — 0, —0.2 and +0.2, respectively. The results for the fit parameter r 
are shown in the table [3]. We have concentrated on A = 00 because this value is the opposite 
extreme to the case A = 0(1/ Np) when ( |1.2| ) follows from the 1/Np expansion. Thus an 



agreement with (1.2) at A = 00 allows to expect that a similar agreement can be found for 



any intermediate value of A. 

The agreement is indeed found for k < 0. At k = (fig. ^ we have r m 0.93. The fits 
are good and although the deviation of r from 1 is larger than in the GN case (fig. [|]), it is 



still similar to that case. The fit in fig. [7| is done for data at k — —0.2 and the results for r 
are the same as in the GN model (within errors). For k = —0.4 the value of r is even nearer 
to 1 than in the GN case (tab. [I], figure not included). There is a clear sign, that /3 is the 
same for k = A = and k < 0, A = oo. 

At k = 0.2, A = oo the fit in fig. |^ shows considerable deviation of r from 1, however. 
The quality of the fit is also poor, i.e. the data cannot be described by the MF Ansatz in 
spite of the presence of the L-dependent parameter r. We suppose that this is an effect of 
the crossover region to the other universality classes at n = k c (\) (discussed in section 6). 

The presented results for the ^-dependence of am F are only a part of the data for amp 
we have accumulated during our studies of the Y 2 models. We have investigated also the 
U(l) symmetric lattice Y 2 model with naive fermions || [K| and for A = and A = 0.5 
performed also simulations of the Z(2) model with various Np. Though less detailed, these 
results are consistent with those presented here. Put together, our results strongly suggest 
that, as y — > at fixed A > and n satisfying — k c (A) < k < k c (A), the fermion mass amp 
behaves according to the scaling law ([1.21). We interpret this as a piece of evidence that the 
Y 2 models belong in the indicated region of A, k parameters to the same universality class 
as the GN2 models with the same Np and chiral symmetry. 

We note that neither the present data nor the MF method of their analysis allow to 
determine the values of the mass gap to compare with the exact result [|15[| . For this purpose 
either the data in the asymptotic scaling region or a much better analysis of finite cut-off 
effects would be needed. 



5 The mass of the XX state 

Needless to say, more observables should be compared if the equivalence of various field 
theoretical models is under discussion. Besides critical exponents, the simplest universal 
quantities are the ratios of masses. A comparison of the spectrum can give valuable hints 
about the universality classes of similar theories. 

Early analytic investigations [O, [TJ| have revealed a rich spectrum in the GN 2 model. 



The mass of the first excited state, which is interpreted as a fermion-antifermion bound 
state, is 

m B = 2m F cos(^- Ty ), (5.1) 

followed by other fermion-antifermion bound states and multifermion bound states. 

We could determine the mass arris in the lattice units. We followed the methods of ref. 
T~9f| and measured the time-dependent correlation function 



Cb{x % ) = ^£M 2 ) , (5.2) 

where x\ is the Euclidean time, x 2 the space coordinate and Mq x the inverted fermion 
matrix. The function Cb{^2) is the correlator (xx{fyxxi x )) a ^ zero spatial momentum. The 
operator xx( x ) is local and therefore very easy to measure, but has small overlap with 
the bound state we are actually looking for. One can write: 

C B (t) = C FF (t) + Z B (er amBt + e - amB{T - t] ) + Z B {-lY (e~ a ™ Bt + e -<^( T -*)) . (5.3) 

Cpp is the free fermion-antifermion cut contribution 

CffH) 4E4 2) (0^)G?(0,4 (5.4) 



U(l) N F =16 2^L/2 
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Figure 9: The ratio mfl/mj? plotted vs. arripL. The results agree for all investigated k,X 
pairs as long as the physical size of the lattice is sufficiently large, arripL > 4. 



G F (0, x) is the free fermion propagator, the corresponding mj? and Zp are determined by 
fitting the measured fermion propagator (x(0)x( x )) ■ 

The mass arris has been obtained at nearly all simulation points. The results for those 
points, at which 

2 < — < £ (5.5) 
amp 2 



are plotted in fig. £| The first inequality in (5T5) selects points without significant lattice 
artifacts, whereas the second one cuts off points with too large finite size effects. In fig. |9] 
we display two physical, dimensionless quantities: the ratio rriB/mp and ampL, the size of 
the system in the units of the fermionic correlation length. 

The finite size effects increase as ampL decreases, and these effects are not expected to 
be universal. But we see in fig. ^ that at least for ampL > 4, all the data are close to the 
expected value ms/Trip = 1.99. This agreement favours the conclusion that for all shown 
(k,X) pairs, including the GN limit, the theory is in the same universality class. 

The large uncertainties introduced by the inexact measurement of arris do not allow to 
draw a definite conclusion, however. The cut contribution Cppit} in eq. ( |5.3|) is large: about 
95%. This makes the determination of not very reliable. We checked, however, that the 
fits for C#(t) with only the cut contribution are significantly worse than those with the full 
expression ( |5.3|) . Nevertheless, we cannot exclude the possibility that the result for the mass 
rriB — Imp is an artefact of the contribution of the fermion-antifermion threshold to Cb(£). 
In any case our results show that the mass of the lightest fermion-antifermion bound state 
is not significantly lower than 2rrip, otherwise the mass arris would be easily measurable. 
This still supports to some extent the universality hypothesis. 
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y 

Figure 10: Two types of continuum limits from the y > in the vicinity of the line k = k c (X). 

6 Transition between the Gross— Neveu and spin model 
universality classes 

We have tentatively concluded that the Y 2 models are equivalent to the GN 2 models in the 
continuum limit y —>■ when k is fixed at some k < K C (X). This still allows that something 
new and interesting happens in the limit when y — > and k — > « C (A) simultaneously. If so, 
this would justify an independent existence of the Y 2 models as field theories with physical 
content different from the GN 2 universality class. 

A hint that this could be the case is provided by the boson spectrum considerations in 
the vicinity of the line k = K C (X) at y = 0. In fig. [H] we indicate the situation at some fixed 
A. The surface y = 0, k < K C {X) corresponds to the high temperature phase of the scalar 
field theory which describes massive scalar particles of mass m^. The continuum limit of 
this theory is obtained (A)_, where am^ —>■ 0. The universality class for A > is 

that of the 2D spin models with the same global symmetry as the scalar field theory. The 
points k < k c (X) are regular points of the scalar field theory, i.e. am^ > 0. Thus the end 
point of the line GN 2 is a regular point and when the continuum limit amp — > is taken 
along this line, the ratio m^/mp becomes infinite. The scalar theory "decouples" and the 
continuum GN 2 model results. 



The situation is different along certain lines of the type X shown in fig. [10]. One can 
imagine a fine-tuned approach towards the line k = k c (A) at y = along which 

TTLfk 

amsi amp — > 0, > a , 

m F 

with some arbitrary constant a. Assuming applicability of the MF approximation, which 
can hold at least for large a and not very close to the point k = k c , one can even calculate 
the corresponding paths (see Appendix). 

Thus the existence of a family of continuum field theories is indicated, labeled by the ratio 
a = m^jmp. Its spectrum might consist roughly of the same states as the GN 2 spectrum 
plus the spectrum of the scalar field theory in the high temperature phase. We see no reason 
why these two "sectors" would not interact with each other, so that these field theories might 
be of a new kind, unexpected when the original action ( |1.1|) of the Y 2 models is considered 
perturbatively. The situation is similar to that found in [ F2"0| , and we suggest that also in 
the case of the Y 2 models there might be "new universality classes which cannot be seen in 
conventional perturbation theory" . 



The new universality classes arising along the lines X could perhaps be also described by 
other actions, different from such that their physical content is more transparent. A 

possible candidate would be the GN 2 model and an additional scalar field theory in the high 
temperature phase, with both sectors coupled to each other in such a way that the spectrum 
of each sector does not change much. As an example one could imagine the coupling of the 
form: ipF(ip)ip with F((p) vanishing faster than if. 

Unfortunately, we see little chance to investigate these issues numerically. They constitute 
a part of the cross-over between two distinctly different universality classes. But the cross- 
over phenomena are notoriously difficult on finite lattices, as one has to work very close to 
the critical point in order to disentangle various mixed scaling tendencies. The difficulties 
we have encountered at A = oo and k = 0.2 suggest that there we are already too close to 
the crossover to be able to apply the simple MF method at those modest correlation lengths 
we could afford. 

7 Discussion and conclusion 

We presented some numerical evidence that the Y 2 models belong to the same universality 
class as the GN 2 models with the same chiral symmetry and Np even if a strong scalar 
self-interaction prohibits the use of the 1/Np expansion. Of course, the universality, or 
equivalence, of field theoretical models requires the equality of all observables like critical 
exponents, mass ratios, etc., in the continuum limit. Thus our evidence, based mainly on 
the scaling behaviour of the fermion mass, might seem to be only fragmentary. However, the 
significance by our results is amplified by the agreement with the MF method. This method 
makes predictions of the same type as the 1/Np expansion in the lowest order, but has been 
found to be applicable in a much wider range of parameters. The MF method thus seems to 
be the suitable analytic method for studying the Y 2 models, provided the pure scalar sector 
is taken into account fully, which is not difficult numerically [0, |18| . The universality is a 
straightforward prediction of this approach. 

We found that the scaling law (|1.2| ) also holds when the bare kinetic term of the scalar 
field is negative. For k < the MF approximation works even better than for k > 0. We 
think that this is due to a larger distance from the cross-over between the GN 2 and the spin 
model universality classes. As the bare scalar field in 2D is dimensionless, the Y 2 models 
can be chosen with quite arbitrary scalar field self-couplings and the ip 4 term, used in eq. 
( |1 . 1| ) , is only an example. However, we believe that our investigation of the Y 2 models 
with A = oo, i.e. \<f>\ = 1 is generic for a broad class of the Y 2 models with quite arbitrary 
scalar self-couplings. 

All that leads us to the conclusion that the Y 2 models belong to the same universality 
class as the GN 2 models with the same chiral symmetry and Np if ft < ft c - Here k c is the 
boundary of the high temperature phase of the pure scalar model with some rather general 
self-couplings. The GN 2 models seem to be the most economic representatives of these 
universality classes, as they do not contain any ballast of irrelevant parameters. 

However, this conclusion does not imply that the Y 2 models are completely superfluous. 
It could be that the continuum limits taken at k = k c could lead to new universality classes, 
as some our speculations indicate. 



A Appendix 



As mentioned in section 6 one can imagine that the ratio a = m^jmp can be tuned when 
one approaches the critical line k = k c (X), y = 0. The value of a depends on the trajectory 
in the bare parameter space along which the continuum limit is approached. Let us consider 
a trajectory of the type X in figure [H], where k c (X) is approached from below. 

In the U(l) case at A = oo the limit y = is the well-known XY model, with the scalar 
mass vanishing near k = n c and the magnetic susceptibility diverging as: 

X ~ U m T V ~ ex P [K K c - «)""] , v = \ r l ,= t- (A- 1 ) 

The approach to k — k c (\), A = oo at fixed ratio a = m^/mp is described by a function 
y(n). According to the MF prediction fll.2p the fermion mass is given by: 

am F « fi F exp I -— — l — ] . (A.2) 



The ratio a now reads: 



This can be solved for y: 



2p (log(ct/! F ) - logm ) x 



where the dependence of \ an d m</> on k is given by ( A.l ). This formula can be used only as 



a rough approximation for the trajectory. On can expect it to be applicable at most for large 
a and not very close to k = k c (but close enough to apply (|A.1[) ). Under these conditions 
eq. ( |A.4| ) can be simplified by keeping only the leading exponential dependence on (k c — k) 
in the r.h.s.: 

r -b{K c -K)- v . (A.5) 



y 2 (n) ~ -^exp 
log a 
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